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ABSTRACT 

The DebtRank algorithm has been increasingly investigated as a method to estimate the impact of shocks in financial networks, 
as it overcomes the limitations of the traditional default-cascade approaches. Here we formulate a dynamical “microscopic” 
theory of instability for financial networks by iterating balance sheet identities of individual banks and by assuming a simple 
rule for the transfer of shocks from borrowers to lenders. By doing so, we generalise the DebtRank formulation, both providing 
an interpretation of the effective dynamics in terms of basic accounting principles and preventing the underestimation of losses 
on certain network topologies. Depending on the structure of the interbank leverage matrix the dynamics is either stable, in 
which case the asymptotic state can be computed analytically, or unstable, meaning that at least one bank will default. We 
apply this framework to a dataset of the top listed European banks in the period 2008 - 2013. We find that network effects 
can generate an amplification of exogenous shocks of a factor ranging between three (in normal periods) and six (during the 
crisis) when we stress the system with a 0.5% shock on external (i.e. non-interbank) assets for all banks. 


Introduction 

The recent economic downturn has made clear that some substantial features of the present financial markets have not been 
properly considered. Regulators 1-1 and academics 4 pointed out the role played by complexity 5-7 in the little understanding of 
the crisis, and in particular the lack of a quantitative assessment for the level of interconnectedness. It has been increasingly 
recognised that the main and simplest way to quantitatively account for the degree of interconnectedness and complexity of 
financial markets is given by the theoretical framework of complex networks . 5-11 By representing financial institutions as 
vertices of a graph we can identify the systemically important ones with the most central vertices . 12-14 Furthermore, the 
evolution of systemic risk can also be modelled by means of dynamical processes on networks . 15-20 

On the one hand, the use of networks makes the quantification and visualisation of interconnectedness possible; on the 
other hand, and perhaps even more importantly, network effects are also responsible for a more subtle, typically unnoticed but 
crucial effect: the amplification of distress. Indeed, while diversification archived through a higher level of interconnectedness 
reduces the individual risk (in the case of independent shocks), it can however increase systemic risk . 21-25 Nevertheless, there 
is no single topological structure that is the most robust in all situations because market liquidity also matters . 26 All these 
issues are presently considered by regulators 27 and the notion of interconnectedness has already entered the debate on “Global 
Systemically Important Banks” (G-SIBs ). 28 

When the banking system is represented as a network, usually propagation of shocks takes place only with removal of ver¬ 
tices in the system, i.e. only after default events. This is an important mechanism for contagion between counterparties , 14,29-31 
although in practice this channel becomes active only if balance sheets are already quite deteriorated 16 or in combination with 
other contagion channels, such as those due to fire sales and overlapping portfolios. 21,22,32, 33 The DebtRank algorithm 19 was 
introduced precisely to overcome this limitation, and to account for the incremental build-up of distress in the system, even 
before the occurrence of defaults. 

At the “microsocopic” level, every financial institution satisfies a balance sheet identity that links the values of its assets 
and liabilities to a capital buffer, which is meant to absorb losses. Balance sheets of different banks are interconnected and 
therefore the mutual interaction between them is expected to play a major role in the emergence of collective properties, as 
it is usually the case for many diverse complex systems. For example, our result for the stability of the system, i.e. that it 
depends only on structural properties and not on the initial state, is a clear example of a general property that finds applications 
in different domains. 


The original DebtRank 34 helped to shift the attention towards interconnectedness as a crucial driver of systemic risk. ' In 
this paper we show that a similar dynamics can be derived from basic accounting principles and from a simple mechanism 
for the propagation of shocks from borrower banks to lender banks. A limitation of the original DebtRank is that banks pass 
on distress to their creditors only once, leading in some cases to a significant underestimation of the level of distress in the 
system. The dynamics proposed here overcomes this limitation by allowing further propagations of shocks. Perhaps the most 
important point is that we are able to characterise the qualitative behaviour of the system by establishing a crucial link between 
the stability of our dynamics and the largest eigenvalue of the interbank leverage matrix. One of the hallmarks of DebtRank is 
that it allowed regulators to monitor at the same time impact and vulnerability of financial institutions by quantifying in terms 
of monetary value the impact of the received shocks. Hence, we test our algorithm on a dataset of 183 European banks listed 
on the stock market. Our analysis shows that systemic risk has consistently decreased between 2008 and 2013, and that banks 
having the largest impact on the system are also the most vulnerable ones. 


Results 

Model description 

We represent the interbank system as a directed network whose nodes are banks. A link of weight A, 7 from node i to node 
j corresponds to an interbank loan from the lender bank i to the borrower bank j of amount A ,y USD. As such, every node 
is characterised by an internal structure given by its balance sheet (see Methods). On the asset (liability) side we distinguish 
between interbank and external assets (liabilities). The interbank assets of bank i correspond to the total amount of outstanding 
loans to other banks within the system, i.e. Y,jAij, while non-interbank assets are called external assets and denoted by Af. 
For every interbank asset A jj in the balance sheet of bank i there is a corresponding interbank liability L;j = Ajj in the balance 
sheet of bank j. As a consequence, links can be interpreted as connections between specific elements of balance sheets, i.e. 
of nodes internal structure. Each bank i also has external liabilities Lf , which correspond to obligations to entities outside 
the system. The equity £) of bank i is defined through the balance sheet identity as the difference between its total assets 
and liabilities. We say that bank i has defaulted if £)■ < 0, i.e. if its total liabilities exceeds its total assets. This is in fact 
only a proxy for a real default event, which is however a common assumption in the literature on financial contagion (see for 
instance 14,29-31 ). 

We now want to write an equation for the evolution of the equity of all banks which remains consistent with the balance 
sheet identity over time. We first define the set of active banks at time t as the set of banks that have not defaulted up to time t: 


.^(t) = {j : Ej(t) > 0}. (1) 

In the following, we will consider a mark-to-market valuation for interbank assets, while liabilities will keep their face value. 
The idea behind this assumption is that the effect of a bank j being under distress is almost immediately incorporated into the 
value of the interbank assets A, ; - held by a creditor bank i, while the obligations of bank j to bank i do not change. When bank 
j defaults, it defaults on all its interbank liabilities, meaning that its creditors will not recover the money that was lent to j and 
Aij will be zero. As a consequence, the balance sheet identity for bank i at time t reads: 


Ei(t)=Af(t)-L?(t)+ £ AyW-LMO. 

j= l 


( 2 ) 


The reason why the sum involving interbank assets runs over all banks active at time t — 1 is that the information about the 
default of other banks is received by bank i with a delay, and accounted for only at the next time step. 

We next assume a simple mechanism for shock propagation from borrowers to lenders. The idea is that relative changes 
in the equity of borrowers are reflected in equal relative changes of interbank assets of lenders at the next time-step: 


Ay(f +1) 


if 7 e =«/(/-1) 

Aij(t) = 0 


(3) 


where the case j -vJ(t 1) ensures that, once bank j defaults, the corresponding interbank assets A ,y of its creditors will 
remain zero for the rest of the evolution. Suppose, for example, that bank j defaults at time s, i.e. Ej (s ■ 1) > 0, but Ej(s) = 0; 
as a consequence, A, ; (s +1) =0, for all i. At time s + 2, since j £ srf (s), the second case will apply, and A, ; (s + 2) = 0. For 
t > s + 2, obviously, A i; (l) will remain equal to zero. 

By iterating the balance sheet identity (2) and the shock propagation mechanism (3), the contagion dynamics can be 
conveniently cast (see Methods for a detailed derivation) in terms of the relative cumulative loss of equity for bank i: hi(t) = 
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(Ei(0) —Ei(t)) /Ej ( 0 ): 


hi(t +1) = min 


N 


1 ,hi(t) + £ A ij(t) [h j(t) - hj(t - 1)] , 

3 =1 


(4a) 


^ ij (0 


0 


(4b) 


where we call A the interbank leverage matrix. 

The above dynamics resembles the DebtRank algorithm already introduced in the literature. 19 An important difference is 
that in the original DebtRank a bank is allowed to propagate shocks only the first time it receives them. In some cases this 
might lead to a severe underestimation of the losses. Let us suppose that bank i is hit at time t by a small shock, which will be 
propagated resulting in additional small shocks at time t + 1 for its creditors. If the network does not contain any loop bank 
i will not be hit again by any other shock. However, if the network does contain loops bank i might be hit at later times by a 
shock which, depending on how much leveraged its borrowers are, might be far larger than the first one, but it will be unable 
to propagate it. Eq. (4) is more general in the sense that as long as a bank receives shocks it will keep propagating them. In 
fact it can be proved that the two algorithm give the same losses on a certain class of networks (as trees), but, in general, the 
losses computed via the original DebtRank are a lower bound to those computed with (4). More precisely, if we shock a single 
node s, the two algorithms will give the same losses for all nodes r such that a unique path from r to ,v exist. If we shock more 
nodes, the two algorithms will give the same losses for all nodes r such that unique and non-overlapping paths between r and 
all the shocked nodes exist. On all the other cases (4) leads to larger losses (see Methods). 

A crucial feature of the dynamics (4) is that its stability is determined by the properties of the interbank leverage matrix 
A (f). Notably, it is possibile to show (see Methods) that when | /L max |, the modulus of the largest eigenvalue of A (t), is smaller 
than one, the dynamics converges to the fixed point A/;(f) = h(t) — h(t — 1) = 0, meaning that the shock is progressively 
damped in subsequent rounds. In contrast, when |A max | > 1 the initial shock will be amplified and at least one bank will 
default. Remarkably, this happens independently on the properties of the initial shock. After the default, according to (4b), 
A (t) will be modified and the same argument will apply to the new interbank leverage matrix. The dynamics will eventually 
converge when the modulus of the largest eigenvalue of A(f) becomes smaller than one. This explains why, even if the system 
is initially in the unstable phase, the dynamics does not necessarily converge to the state in which all the banks default. When 
a bank defaults it is effectively removed from the system when the interbank leverage matrix is updated. The new, reduced 
system could now be in stable phase, and thus converge to the stable fixed point. The important point here is that, although the 
exact values of final losses will depend on the initial shock, the ability of the system to amplify distress and lead to defaults is 
an exclusive property of the leverage matrix. This result confirms the importance of the leverage matrix for the amplification 
of shocks within the context of systemic stability, as suggested by, 17 albeit for a different contagion mechanism (the so-called 
Furfine algorithm 36 ). 


Application to the European banking system 

We now apply the introduced algorithm to the European banking system. We use data from the balance sheets of 183 publicly 
traded European banks between 2008 and 2013. Available data only contain information about the total amount of interbank 
borrowing and lending for each bank, which are respectively the sum over rows and columns of the matrix of interbank assets 
Aij. Therefore, we resort to a two-steps reconstruction technique 37-39 to infer plausible values for all the entries of the matrix. 
In the first step we build the topology of the network using a so-called fitness model, while in the second step we assign 
weights to links using the RAS algorithm 40 (see Methods for more details about the reconstruction procedure). Due to the 
stochasticity of the first step, we sample 100 different networks, which will be used in the following experiments. 

As a first scenario, we consider a shock affecting all banks simultaneously at time 1 = 1 corresponding to a relative deval¬ 
uation a of their external assets. Following 17 we measure the response of each bank to the shock in terms of its contribution 
Hj(t) to the relative equity loss of the system: 


Hi(t) = 


Ej(0) 
LiE,( 0) 


hi(t) 


Ei( 0) 

LiEi( 0 ) ' 


(5) 


The direct effects of the shock in terms of relative equity loss are //,• (1), while the effects of contagion are computed using 
the algorithm introduced here, which is run until convergence (see Methods for more details). In the top panels of Fig. 1 
we compare the total relative equity loss H(t ) = £,//,(/) directly due to the initial shock (i.e. at time t = 1) with the one 
that includes losses generated by the contagion dynamics (i.e. at the convergence of the algorithm), for all the years, and for 
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Figure 1 . Relative equity loss for the system of 183 publicly traded EU banks between 2008 and 2013. All banks are 
subject to an initial shock consisting in the devaluation of their external assets by a factor a. The violet curves represent the 
relative equity loss that is directly due to the initial shock, while the blue curves include losses due to the contagion dynamics. 
Every point is the average over 100 reconstructed networks and the semi-transparent region covers the range between the 
minimum and maximum across the sample, a is fixed in the top panels and equals to 0.5% (panel a) or 1% (panel b). We see 
that the amplification effect is reduced from 2008 to 2013. Bottom panels refer to 2008 (panel c) and 2013 (panel d). We see 
that the relative equity loss saturates for large enough shocks. In 2008 the saturation already occurs for shocks as large as 


0.5%. 


a = 0.5% and 1%. The overall behaviour resembles the one reported in 37 obtained using the original DebtRank. However, as 
already discussed, the relative equity losses observed here are larger by a factor ranging from 1.3 in 2008 to 1.7 in 2013. 

We further test this scenario in the bottom panels of Fig. 1 by focusing on 2008 and 2013 and letting a vary between 
0.5% and 5.5%. The relative equity loss experienced by the system increases as we increase a, until it reaches a saturation 
point. For large enough values of a most of the equity of banks is already wiped out by the initial shock, implying that the 
amplification due to the contagion dynamics decreases with a. Interestingly, we observe that in 2008 the amplification pushes 
relative losses of equity to saturation levels already for values of a as small as 0.5%, while in 2013 shock five times larger are 
needed to reach similar relative losses. 

As a second scenario, we consider the case in which a single bank at a time is shocked, a shock still being a devaluation of 
its external assets by a relative amount a, and the experiment is repeated for each bank. The idea is to decompose the systemic 
importance of a bank into its impact on the system and into its vulnerability with respect to shocks affecting other banks. We 
then proceed to define the impact of bank i as the relative equity loss of the system when bank i is shocked. Instead we take 
as a measure of its vulnerability the average of lu(t ) over all the experiments. We then rank banks in descending order both in 
terms of impact and in terms of vulnerability and present the results for 2008 and 2013 and for a = 0.5% in form of a scatter 
plot in Fig. 2. We can see that the most dangerous banks, i.e. the banks having the largest impact on the system, are also the 
most vulnerable ones. 
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Figure 2. Scatter plot of impact and vulnerability (reverse) rankings in 2008 (panel a) and 2013 (panel b). An initial 
shock corresponding to a 0.5% devaluation of its assets is applied to one bank at a time, and the experiment is repeated for 
each bank. The impact of a bank is measured as the relative equity loss experienced by the system when that bank is shocked. 
The vulnerability of a bank is its relative equity loss averaged over all the experiments. In addition, we average impact and 
vulnerability across a sample of 100 reconstructed networks. Finally, we build reverse ranking (i.e. in descending order) of 
both quantities, so that larger values on both axes correspond to more impactful and more vulnerable banks. Bubble size is 
proportional to the total assets of the corresponding bank. The most dangerous banks are also the most vulnerable. 


Discussion 

By iterating the balance sheet identity we derive an equation for the evolution of banks’ equities. We then consider a general 
shock propagation mechanism in interbank networks so that the value of interbank assets of lenders depends on the level 
of distress of their creditors. The resulting dynamics is closely connected with the DebtRank algorithm recently introduced 
in the literature 19 as an effective shock propagation dynamics and it provides a clear economic intuition for its dynamical 
variables, in terms of basic accounting principles. We prove that, in general, the original DebtRank gives a lower-bound for 
losses computed with our methodology, but, for a certain class of shocks, the two algorithms are equivalent on trees. More 
importantly, we show that the capability of the system to amplify an initial shock depends only on the modulus of largest 
eigenvalue Amax of the matrix of interbank leverages: When A max < 1, additional losses induced by subsequent rounds of the 
dynamics are attenuated over time. In contrast, when A max > 1 a small shock will be amplified and cause at least one bank 
to default. This finding can be important from a regulatory perspective, as one could monitor the evolution of Amax over time 
to check if the system is entering the unstable regime. 

To showcase our algorithm, we apply it to a system composed of 183 European publicly traded banks. We characterise 
the response of the network to different shock scenarios. Our analysis shows that the amplification of shocks due to interbank 
contagion consistently decreases from 2008 to 2013, and that in 2008 small shocks are enough for all banks to be significantly 
distressed. By performing stress tests in which banks are initially shocked one at a time, we are able to compute both the 
impact of a single bank on the system and its vulnerability to shocks initiated by other banks. From a systemic standpoint, 
it would be desirable that systemic impact and vulnerability were anti-correlated, so that the most dangerous bank are also 
the most robust, and vice versa. In fact, this does not happen: our analysis shows that the most dangerous banks are also the 
most vulnerable, meaning that systemic risk is concentrated in a few key players, which should therefore be the objective of 
effective macroprudential regulation policies. 


Methods 

Balance sheet basics 

A balance sheet summarises the financial position of a bank. It consists of assets, which have a positive economic value 
(e.g. stocks, bonds, cash), and liabilities, which are obligations to creditors (e.g. customers’ deposits, other debits). The 
difference between the value of assets and liabilities is called equity, and the following (balance sheet) identity holds: assets = 
equity + liabilities. A bank is said to be solvent as long as its equity is positive. Once a bank is insolvent, even if it sold the 
entirety of its assets, it would not be able to repay its debts. As a consequence, we use insolvency as a proxy for default. 
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Model dynamics 

The equation for the evolution of the cumulative relative loss of equity hi(t) = (£)(0) — Ej(t)) /E;(0) can be derived from the 
balance sheet identity. From (2), supposing that (i) external assets and liabilities do not change, (ii) interbank liabilities are at 
face value and also do not change, and (iii) interbank assets are marked-to-market, 

E i {t + \)-E i {t)= £ Aijit +1) — £ Ajjit) 

= E [Ay(f + 1)-Ay(0]- E Aij[t+ 1), 


where in the second line we have isolated a potential contribution coming from the nodes that were active at time t — 1, but 
became inactive at time t. Using (3), we see that the last term in (6) vanishes, so that we have: 


E i (t+l)-E i (t)= E 


Aijit) 


jes/(t~\) Eji 1 !) 
Ay(0) 


[Ej(t)-Ej(t-l)} 


(7) 


je J( t . x) E j^) 


[Ej(t)-Ej(t- 1)] 


where in the second line we have recursively applied (3) and used Ay(l) = Ay(0) (only equities change at time t — 1, assets 
start to change at time t = 2). We can now define the matrix A: 


Aijit) 


T0T ifj € #?(t-l) 
0 if j srf{t — 1) 


( 8 ) 


and write the equation for the evolution of equity: 


Ei(t + 1) = max 


N 


o,£,-(0 + E^-(0[^-(0-^-1)] . 

j =1 


(9) 


where the max accounts for the fact that once a bank defaults its equity cannot become negative. From (9), it easily follows 
that: 


hi(t + 1) = min 


N 


Mi(0 + E Aijit) [hj{t ) - hj(t - 1)] 

;=i 


( 10 ) 


where Ajjit) = Ajj{t)Ej(0) /£,-(0), so that A(/) can be interpreted as a reduced interbank leverage matrix, where columns 
corresponding to banks defaulted up to time t — 1 have been set equal to zero. As the equity of defaulted banks does not 
change anymore after reaching zero, the rows of the leverage matrix corresponding to defaulted banks can be set equal to zero 
too. 


Relation to DebtRank 

The original DebtRank 19 has the following dynamics: 


hjit + 1) = min 


1 ,h{{t)+ E Wjjhjit) 


= min 


1 ,hi(t)+ E Wij[hj(t)-hj{t-l)] , 


( 11 ) 


where Wy = min(l,Ay), and srf'{t) = {j : hjit) > 0 and hj(t — 1) = 0}, and the last term in the second line can be added 
because it is always equal to zero. Let us note that the definition of £f'{t) implies a different stopping criterion. In fact, in the 
original DebtRank nodes propagate shocks only once, immediately after the shock has been received. In our setting, instead, 
they could propagate shocks until they default. There are two main differences with respect to (4a): (i) the summation in (11) 
involves less terms than the summation in (4a) since srf'it) C &/(t) C g/(t — 1) (the set of active nodes becomes smaller and 
smaller as banks default); (ii) Wy < Ajj, for all i and j. As a consequence, (11) provides a lower bound to relative cumulative 
losses of equity computed with (4a). 
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In order to understand the role of the network topology, let us focus our attention on a node r. From (3) we see that a 
shock can reach r only through the neighbours r borrows from, which in turn can be reached by a shock only through the 
neighbours they borrow from. In other words, if a single node s is shocked at some time t, the only possible way for r to 
experience the effects of such shock (at later times) is that a path from r to s exists. Let us for a moment suppose that such 
path r — > i\ —> 4 ■ ■ ■ — > ip -1 — > s is unique (and of length p)\ then there will be also a unique path leading from any node 4 
to s. The shock will propagate to node i p \ at the time t + 1, but, if no additional node is shocked, and since no additional 
paths exist between s and i p -\, the status of node i p \ will not change from time t + 1 to time 1 + 2. Similarly the status 
of node z„_ 2 will change only at time 1+2, and so on, until the shock reaches node r at time t + p. The status of any node 
4 on the path will change only at one time step. As a consequence, the result will be the same as if each node were active 
only when reached for the first time by the shock, as in the original DebtRank. However, this is true only for the nodes r 
such that a unique path connecting them to the only shocked node s exists. If there are additional paths between r and .v the 
shock will propagate also along those paths, resulting in additional losses at the node r. In particular this is trivially true if the 
subgraph of nodes reachable (backwards) from the shocked node s is a tree. If more than a single node is shocked, and if r is 
reachable (backwards) from more than one of them, then, even if the graph is a tree, the (cumulative) loss experienced by r at 
the end could be larger than if a stopping criterion a la DebtRank were used. In particular the loss will be larger if the paths 
are overlapping, while it will be equal if the paths are not overlapping. 

Stability properties 

Let us assume for simplicity that no banks default during the whole evolution, so that A is constant over time (see (8)). 
Defining A h(t) = h(t) — h(t — 1), (4a) can be written in matrix notation: 

Ah(t + 1) = AAh(t) 

= A‘Ah{l)=A‘h{l), (12) 

as A/z(1) = /z(l) — /z(0), and /z(0) = 0. By summing over all the time steps up to t + 1 one gets: 

t +1 f+i 

/z(f + l) = £A/z(s) = £A s /z(l). (13) 

5=0 5=0 

A/z = 0 is always a fixed point of the map (12), and it is stable as long as the modulus of the largest eigenvalue A max of A is 
smaller than one, meaning that the dynamics will damp subsequent propagations of an initial shock over time. In this case the 
sum in (13) will asymptotically converge to: 

/z°° = (1 — A) -1 /z(l). (14) 

In contrast, if | A max | > 1, A/z (t ) will become increasingly larger, leading to the default of at least one bank, independently from 
the initial shock. 

Eq. (12) clearly describes the first stages of the dynamics, up to the first default. Nevertheless, since the reduced leverage 
matrix does not change between two subsequent defaults, ( 12) also holds between one default and the next one, provided that 
A is replaced with the correct reduced leverage matrix A(t). As a consequence, the dynamics will remain explosive as long 
as the modulus of A max (f), the largest eigenvalue of A(f) is larger than one. As more and more banks default lAmaxCOI will 
eventually become smaller than one, and the dynamics will finally converge. 

It should be noted that modifying the original DebtRank dynamics (11) by allowing banks to propagate shocks as long as 
their equity is positive would lead to a double-counting of losses. Let us suppose again for simplicity that no banks default 
and that W = A. Iterating (11) leads to h(t + 1) = (1 +A) r /z(l), and this quantity is always larger than the one obtained from 
equation (13), i.e. L x . = qA s /z(l). 

Data 

For our analysis we use the same dataset used in. 37 Information on banks’ balance sheets are taken from the Bureau Van Dijk 
Bankscope database for 183 European banks that were publicly traded between 2008 and 2013. From this data source we 
extract information about: equity, total assets, total liabilities, total interbank assets A,- and total interbank liabilities Lj. For 
details about the handling of missing data, the reader should refer to the aforementioned reference. ’ 7 

As mentioned in the main text, the procedure to reconstruct a matrix of interbank assets A, ; develops in two steps: in the 
first we generate a binary adjacency matrix, which encodes the topology of the network. This is done via a fitness model, 41 

^out^in 

conveniently modified for directed networks. A link from bank i to bank j is inserted with probability p,j = | ' 1 , IUI ; n . , where 

l+ZX; X ■ 
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the fitness values of each banks are computed as x° ut = A ,7 Aj and x'" = L,/ ; Lj , and the parameter z is fixed to attain the 
desired network density (the number of links in the network divided by the number of possible links). In this paper we have 
set z so that the density of the network is 5%. We then draw 100 networks according to the probabilities pij. For each network 
thus obtained, we then proceed to assign weights Ag to the links. To this end we use the RAS algorithm. 40 This consists in 
the iteration of a map whose n-th step is: 


A (") — 

A ij ~ 


4 "* 

iaST 


-A, 


• (n+ 1) 


L-A 


(«) 


Li. 


At convergence, the above iteration ensures that / ,4 J; = A, and A (/ = L j for all banks. Obviously, one must have that 
A; = Y., L,. Since this is not the case for our data, we rescale liabilities Li so that the above relation holds. 
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